Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 19 November 2010 



(MN style file v2.2) 



O 



Detection of a radial velocity gradient in the extended 
local disc with RAVE 

o": 

O - A. Siebert 1 *, B. Famaey 1 , I. Minchev 1 , G.M. Seabroke 2 , J. Binney 3 , B. Burnett 3 , 

K.C. Freeman 4 , M. Williams 5 , O. Bienayme 1 , J. Bland-Hawthorn 6 , R. Campbell 5 , 
§ J.P. Fulbright 7 , B.K. Gibson 8 , G. Gilmore 9 , E.K. Grebel 10 , A. Helmi 11 , U. Munari 12 , 
£ : J.F. Navarro 13 , QA. Parker 14 , WA. Reid 14 , A. Siviero 5 ' 12 , M. Steinmetz 5 , F. Watson 15 , 
R.F.G. Wyse 7 , T. Zwitter 16 - 17 

1 Observatoire Astronomique, Universite de Strasbourg, CNRS, Strasbourg, France 

2 Mullard Space Science Laboratory, University College London, Holmbury St Mary, Dorking, RH5 6NT, UK 

■ 3 Rudolf Peierls Centre for Theoretical Physics, Oxford, UK 
E 4 Australian National University, Canberra, Australia 

| 5 Astrophysikalishes Institut Potsdam , Potsdam, Germany 
• , 6 Sydney Institute for Astronomy, University of Sydney, Sydney, Australia 

7 Johns Hopkins University, Baltimore, MD, USA 

8 University of Central Lancashire, Preston, UK 

9 Institute of Astronomy, Cambridge, UK 

10 Astronomisches Rechen- Institut, Zentrum fur Astronomie der Universitat Heidelberg, Heidelberg, Germany 

11 Kapteyn Astronomical Institut, Groningen, the Netherlands 

12 Astronomical Observatory of Padova in Asiago, Asiago, Italy 

13 University of Victoria, Victoria, Canada 

14 Macquary University, Sydney, Australia 

r- H ■ 15 Anglo- Australian Observatory, Sydney, Australia 

^ ■ 16 University of Ljubljana, Faculty of Mathematics and Physics, Ljubljana, Slovenia 
£\\ ' 17 Center of excellence, SPACE-SI, Ljubljana, Slovenia 

C\ ' 
O ; 

\l Accepted . Received ; in original form 19 November 2010 

ABSTRACT 

Using a sample of 213 713 stars from the Radial Velocity Experiment (RAVE) sur- 
vey, limited to a distance of 2 kpc from the Sun and to \z\ < 1 kpc, we report the 
detection of a velocity gradient of disc stars in the fourth quadrant, directed radially 

■ from the Galactic centre. In the direction of the Galactic centre, we apply a simple 

method independent of stellar proper motions and of Galactic parameters to assess 
the existence of this gradient in the RAVE data. This velocity gradient corresponds 
to \K + C| >3 kms _1 kpc _1 , where K and C are the Oort constants measuring the 
local divergence and radial shear of the velocity field, respectively. In order to illus- 
trate the effect, assuming a zero radial velocity of the Local Standard of Rest we then 
reconstruct the two-dimensional Galactocentric velocity maps using two different sets 
of proper motions and photometric distances based either on isochrone fitting or on 
iv"-band magnitudes, and considering two sets of values for the Galactocentric radius 
of the Sun and local circular speed. Further observational confirmation of our finding 
with line-of-sight velocities of stars at low latitudes, together with further modelling, 
should help constrain the non-axisymmetric components of the Galactic potential, 
including the bar, the spiral arms and possibly the ellipticity of the dark halo. 

Key words: Stars: kinematics - Galaxy: fundamental parameters - Galaxy: kine- 
matics and dynamics. 



1 INTRODUCTION 
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Our Milky Way Galaxy is a unique laboratory in which 
to study galactic structure and evolution from stellar kine- 
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matics. The story of the efforts to obtain such kinematic 
data nicely illustrates how theoretical progress and data 
acquisition have to go hand in hand if one wants to gain 
insight into the structure and history of the Galaxy. Un- 
til recently, most studies have been limited to the so- 
lar neighbourhood. The zeroth order approximation for 
modelling these kinematic data assumes an axisymmet- 
ric model, in which the Local Standard of Rest (LSR) 
is on a perfectly circular orbit. However, the local veloc- 
ity field in the solar neighbourhood already displays pos- 
sible signatures of the non-axisymmetry of the Galactic 
potential in the form of stellar moving groups contain- 
i ng stars of ver y different ages and c h emical compositions 
dDehnen I [l998l: IFamaev et al. I 120051: iQuillen fe Minchev I 



120051 : iBensbv et al. 1120071 : iMinchev et al~ Il2010h . 



In external galaxies, large-scale distortions of the (line- 
of-sight) isovelocity contours of the gas have l ong been recog - 
nised in two-dimensional velocity maps (e.g. iBosma II1978T ). 
allowing a quantification of the strength of non-circular 
strea ming motions (e.g., ISellwood &; Zanmar Sanchez I 
2010 ). A prominent exam ple is the galaxy M81 



(e, 



Adler Sz W efstpfahl 1996), which has roughly the same 



circular velocity as the Milky Way and in which the radial 
streaming motions of the gas reach maxima of the order 
of 50 km s _1 . However, the stellar kinematics often turn 
out to be slightly m ore regular and sym metric than the 
gas kinematics ( e.g. IPizzella et al~l 120081 ). For instance, 
iRix &; Zaritskv I (| 19951 ) estimate from 18 face-on spiral 
galaxies that the azimuthally averaged radial streaming 
motions of disk stars should typically amount to a maxi- 
mum of 7 kins" 1 due to lopsided distortions and 6 kins" 1 
due to intrinsic ellipticity. 



In the Milky Way, non-axisymmetric motions are 
clearly visible in radio-frequency observations of gas flow 
in the inner Galaxy, implying that the amplitude of spiral 
structure could be larger by a factor ~1.5 tha n its amplitude 
in the near-infrared luminosi ty density (|Bissantz et al. I 
120031 : IFamaev fe Binnev I l20o3 ) . On the other hand, radio 
interferometry used to determine the parallax of star- 
forming regions in the Perseus spiral arm, together with 
the measurement of their line-of-sight (los) velocities, show 
that these star-forming regions have peculiar non-circular 
motions amounting to ~20 km s _1 , assuming standard 
IAU Galactic parameters and a nearly-flat rota tion curve 
(|Xu et al. 112005 iBinnev lliooS iReid et al. ll2Q09h . However, 
the amplitude of such non-circular mot ions strongly depends 
on th e Galactic parameters adopted ([McMillan &; Binnev I 
l20ld ). namely the Galactocentric radius of the Sun Ro, 
the local circular speed v c o, and the reflex solar motion, as 
well as the model for the rotation curve. Another example 
of such non-circular motion is the 3 kpc arm which has 
^ios = — 80km s _1 towards the Galactic centre (GC) and 
is known to have a counte rpart on the far side of the GC 
(|Dame &; Thaddeus 1 120081 ). Its origin is associated to the 
Galactic bar showing that non-circular motion of the similar 
amplitude can be caused by various types of perturbations 
and not necessarily by spiral arms. 



Here, we use line-of-sight velocit}Q data (v\ os ) from the 
RAVE survey to investigate whether the mean radial motion 
of disc stars remains zero in the extended solar neighbour- 
hood as compared to the one in the solar neighbourhood, 
or whether large-scale radial streaming motions are present 
in the vicinity of the Sun, which would translate in a ra- 
dial velocity gradient in one or the other direction. In Sec. 2 
we present the sample used for our analysis while in Sec. 3 
we present the detection of a radial velocity gradient using 
RAVE 

^ios 5 as well as the associated Oort constants and 2- 
dimensional velocity field. Finally in Sec. 4 we discuss the 
implications of our results and present our conclusions. 



2 THE RAVE SAMPLE 

Our analysis is based on the RAVE catalogue 
(|Steinmetz et al. lliooi : IZwitter et al. 1120081 ) which provides 
v\os to 2 kms" 1 and stellar atmospheric parameters and 
distances to 30% for a large number of bright stars in the 
southern hemisphere with 9 < I < 12 using moderate 
resolution spectra (R^7500). RAVE selects its targets 
randomly in the /-band interval, and so its properties are 
similar to a magnitude limited survey. The RAVE cata- 
logue is cross-matched with astrometric (PPMX, UCAC2, 
Tycho-2) and photometric catalogues (2MASS, DENIS) to 
provide additional proper motions and magnitudes making 
RAVE a suitable tool to study the M ilky Way. 

We use the latest version of the I Z witter et al. I (|2010h 
catalogue, which contains four sets of distances in addition 
to the standard RAVE data (line-of-sight velocities and stel- 
lar parameters). These four sets of distances are (i) the dis- 
tance moduli computed using three sets of isochrones (in- 
cluding the Pa d ova is ochrones), and (ii) the distances of 
iBreddels et"al~l (|2010l ). In this paper we will use the dis- 
tance moduli based on the Padova isochrones but our re- 
sults do not strongly depend on the set of isochrones used 
(see also Sect l3.2l where all distances are multiplied by an 
arbitrary factor). This catalogue is cleaned from potentially 
problematic fields using a list kindly provided by G. Matije- 
vic (private communication) and we further select stars with 
\z\ < 1 kpc to minimise contamination by halo stars. This 
leaves us with a sample of 213 713 stars with known dis- 
tances, atmospheric parameters (logg, T e ff, [m/H]), proper 
motions as well as v\ os measurements that can be used to 
compute the 3D space velocities. 

To verify that our analysis does not depend on the 
method used to determine the distances or on the source 
of proper motions, we also use stars for which distances 
and proper motions are obtained independently. To do so, 
we select a sub-sample of red clump candidates for which 
the distance can be estimated using the if -magnitude alone 
because of the narrow, almost Gaussian, luminosity func- 
tion of the red clump. Red clump candidates are selected 
in the (logg, J — K) plane, as illustrated in Fig. [1] We 
select the over-density using 0.5 < J — K < 0.735 and 



1 Throughout this paper we will use the notation v\ os or line-of- 
sight velocity to refer to the heliocentric radial velocity measured 
by RAVE while radial velocity or Vr means galactocentric radial 
velocity 
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Figure 1. Distribution of stars in the (log g, J—K) plane near the 
red clump for the RAVE internal release. Colour coding follows 
the number of stars per bin of 0.03 mag width in J — K and 
0.05 dex in logg. The distribution is smoothed over 3 pixels to 
reduce the noise. Stars in the white box are identified as red-clump 
stars. 

1.8 < logg < 2.8: the red limit helps reduce the contam- 
ination by the upper red-giant-branch stars. Distances are 
obtained assuming all stars in the sample are red clump 
stars using a Gaussian luminosity function for this popula- 
tion with Mr = —1.6 and ok — 0.22. We apply the same cut 
in z as for the full sample. This results in a sample containing 
29 623 red clump candidates. More details o n the contam- 
inatio n and distance errors can be found in ISiebert et al. I 
(2008). We cross-match this sample with the UCAC3 cat- 
alogue to provide a different source for the proper motion 
- the RAVE catalogue contains predominantly Tycho-2 and 
PPMX proper motions. Hence, for this sample, the distances 
and proper motions are independent from the full RAVE 
sample. The mean height above/below the Galactic plane 
for each sample is shown in Fig. [21 



3 A RADIAL VELOCITY GRADIENT 

3.1 Line-of-sight velocity detection 

Most of the RAVE stars lie in the fourth Galactic quadrant 
with the tails of the distribution in Galactic longitude ex- 
tending in the first and third quadrants. For our analysis to 
remain independent of Galactic parameters and proper mo- 
tions, we first look in the GC (and anticentre) direction with 
|/| < 5° (and 175° < I < 185°). In this case, the Galacto- 
centric radial velocity is close to the measured line-of-sight 
velocity projected on the plane (or its opposite in the an- 
ticentre direction), with only a small contribution from the 
proper motion vector. 

In Fig. [3] we plot the projection onto the plane of the 
mean v\ os in bins 200 pc wide as a function of d cos £ cos b. 
The left panel is for the full sample and the right panel 
for the red clump candidates. We then compare the ob- 
served mean velocities to the expected velocities for a thin 
disc in circular rotation at v c o and adding a radial gradi- 
ent d{VR)/dR =3, 5, and 10 km s _1 kpc -1 , where Vr is 
the Galactocentric radial velocity (see full curves from top 



to bottom in each panel in Fig. [3]). The curves are drawn 
using 5 000 Monte Carlo realisations of the samples, replac- 
ing the measured velocities by random realisations of a thin 
disc ellipsoid. The model used is the same for both sam- 
ples, kinem atic parameters for the thin disc ellipsoid being 
taken from iBinnev fe Merrifieldl (|l998T ) table 10.4. Clearly, 
our data are not compatible with a disc in circular rotation 
because the mean projected line-of-sight velocity is system- 
atically lower than the expected mean velocity. The offset 
between the two panels reflects the different distribution in 
R of the two samples. 

Because our samples are not perfectly oriented towards 
the GC, the rotation component could bias our measure- 
ment. Therefore, we reproduce this by modelling a thick disc 
population which lags the circular rotation by 36 km s _1 . 
The result for a thick disc is presented as dashed curves 
and is computed only in the case of circular rotation. The 
results do not differ significantly from the thin disc results 
indicating that the contribution from v c o to our line-of-sight 
velocities is marginal and hence that our result truly mea- 
sures a radial velocity gradient in the RAVE data. 



3.2 Oort constants 

As first worked out by Oort (1927) for an axisymmet- 
ric Galaxy, and later generalised to non-axisymmetry (e.g., 
Chandrasekhar 1942), a practical way to study the global 
streaming motion of field stars in the Galaxy is to con- 
sider an arbitrary smooth velocity field and Taylor expand it 
around the position of the Sun. Such a Taylor expansion of 
the Cartesian velocity components (U,V) to first order yields 
the four Oort constants A, B, C, and K, measuring the az- 
imuthal (A) and radial (C) shear, the local divergence (K), 
and the vorticity (B) of the velocity field (see e.g. Kuijken & 
Tremaine 1994, Oiling & Dehnen 2003, hereafter KT94 and 
OD03 respectively). 

While we note that such a first order Taylor expan- 
sion might not be a valid description of the velocity field 
for sample depths larger than 1 kpc (as already noted in 
Pont et al. 1994 and OD03), we also note that, at least in 
the direction of the GC, a constant linear velocity gradient 
OVr/OR = K + C between —3 and —10 km s _1 kpc -1 is suf- 
ficient to reproduce the trend seen in Fig|3] In view of this, 
we compute the Oort constants in the classical way using 

( Vrad / r ) = K + A sin 2£ + C cos 21 , 

(vt/r) = B + A cos 2£ - C sin 2£ , (1) 

where v ra d (resp. v t ) is the heliocentric radial (resp. 
transverse) component of the velocity vector in the Galactic 
plane corrected for the motion of the Sun with respect to 
the LSR. r is the distance from the star to the Sun in the 
Galactic plane and ABC and K are the Oort constants 
we try to measure (see also KT94, OD03). We then fit the 
right-hand side of eq. [1] to the left-hand side that is built 
upon our observations with v ra d depending on v\ os and fib 
while v t depends on ji£. 

We restrict the sample to the —140 < £ < lOdeg in- 
terval, the noise in the computed mean velocities prevent- 
ing a reliable fit outside of this interval. In order not to 
rely too much on the proper motions and mostly on RAVE 
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Full sample 



Red clump candidates 





Figure 2. Mean \z\ distance to the Galactic plane in kpc for the two samples studied: full sample (left), red clump sample (right). The 
two panels contain 60 x 60 pixels in a box 4x4 kpc in extent. Both samples have an upper limit on \z\, max(\z\) = 1 kpc. The Sun is 
located at the centre of the map, X is increasing in the direction of the Galactic centre and Y towards the Galactic rotation. The solar 
circle is shown as a dashed line. The circle in the centre represents a sphere 125 pc in radius. 
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Figure 3. Projection of the mean RAVE -ui os on the Galactic plane in distance intervals of 200 pc in the direction of the Galactic centre 
(|/| < 5°) and anticentre (175° < I < 185°) for our two samples. Left: full sample containing 10 682 stars. Right: red clump sample 
containing 1 328 stars. The open circles are our data with 1 — a error bars. The dashed curve corresponds to a thick disc model with no 
net radial motion, while the full curves are for a thin disc with a radial velocity gradient of 0, 3, 5 and 10 km s — 1 kpc -1 from top to 
bottom. 



data, the median proper motion error of the sample being 
1.8 mas yr _1 , we first fit only the constants A, C and K 
to the line-of-sight velocity in the Galactic plane v ra d, thus 
using only v\ os and fib • We then use all the available in- 
formation, and thus fit also the constants A, B and C to 
the tangential velocity in the plane (based on in). Once the 
Oort constants are fitted to the mean velocity field of the 
observed population, we correct them for the asymmetric 
drift as in OD03: since the RAVE catalogue is dominated 
by the old thin disk, we use an asymmetric drift of 20 km 
s _1 . The results are listed in Table 1, and compared to those 
of KT94 and OD03. We note that, although we used the up- 

I 1 A k 

dated values of Sc honrich et al. I ([2010 ) for the reflex motion 
of the Sun, our value for C, especially when fitting only the 
line-of-sight velocity in the Galactic plane, is in remarkable 
agreement with the value found by OD03 (after the mode- 



mixing correction) from proper motions of old red giantqj 
in the ACT/Tycho-2 catalogues. However, our value of C 
(and K) is in disagreement with the one reported in KT94 
from a review of the older literature. Let us finally note that 
the formal errors quoted in Table 1 were computed using a 
Monte- Carlo sampling but are large underestimations of the 
true errors, due to possible systematics including those listed 
in OD03, as well as the absence of a full longitude coverage. 
This underestimation of the errors is confirmed by the fact 
that, when taking or not taking into account the longitudi- 
nal proper motions, the values of the Oort constants A and 
C differ by more than 6 formal sigmas and 4 formal sigmas, 



2 OD03 found different values for different populations (espe- 
cially, lower \C\ for lower velocity dispersion populations), but 
their best tracers of the "true" Oort constants are the red giants, 
also affected by an asymmetric drift of roughly 20 km s _1 
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respectively. The fit to v ra d alone (3 rd line of Table 1) is 
more reliable as it mostly relies on the RAVE line-of-sight 
velocities and is of better quality than the one including lon- 
gitudinal proper motions (the reduced % 2 values are 1.4 and 
2.6 respectively). 

This fitting procedure of the Oort con- 
stants thus confirmed our finding of Sect 13. II that 
\dV R /dR\>3 kms _1 kpc _1 in the RAVE catalogue. 
Let us note that, as we rely on a priori known distances, our 
result is not affected by mode-mixing. However, in order to 
check whether the qualitative conclusion of the existence 
of non-zero velocity gradient depends on our estimates of 
the distances, we assumed that all distances were wrong 
by a factor /, and re-applied the same procedure. Results 
for / = 0.8 and / = 1.2 are listed at the end of Table 1, 
confirming that K + C ^ in both cases. 

3.3 Two-dimensional velocity field 

To further investigate the velocity gradient detected with 
RAVE, we now compute the full 2D velocity field for our 
samples (making use of both the line-of-sight velocities and 
the proper motion vectors). 

The velocity fields in the Galactocentric reference frame 
are computed using 

Vgsr = V* /0 + Vq/lsr + Vlsr/gsr, (2) 

where Vq/lsr is the Sun's velocity with respect to the lo- 
cal standard of rest (LSR) and V lsr/gsr is the motion 
of the LSR with respect to the GC. Obviously, a possible 
radial motion of the LSR itself can be a priori difficult to 
disentangle from large-scale streaming motions. In order to 
circumvent this, we use hereafter the standard assumption 
that the LSR is on a circular orbit0. The velocity of stars 
Vgsr is hereafter decomposed into three components in the 
axisymmetric galactocentric reference frame (Vr, V^, V z ), V z 
pointing towards the south Galactic pole. 

This projection depends on two parameters : Ro, 
the distance to the centre of the Galaxy and v c o, the 
local circular speed . Unf ortunately, as pointed out by 
iMcMillan fc Binnev I <|2010h , neither of these parameters are 
well constrained. The ratio v c o/Ro is however better con- 
strained, lying in the range 29.9 — 31.6 km s _1 kpc -1 with 
some im pact on the measured V_qjl sr as shown in table 1 
and 2 of McMillan fe Binnev I <|2010h . 

Therefore, we compute the 2D velocity maps for two 
se ts of parameters. In both cases we fix V /lsr to the values 
of ISchonrich et al~l (|2010h . The first set of parameters uses 
the standard IAU values v c o — 220 km s _1 and Ro — 8 kpc 
while for the second s et we use the model 2 from table 2 of 
IMcMillan fe Binnev I (|201Ch . v c0 = 247 km s" 1 and R = 
7.8 kpc, which h as a solar motion vector close to the latest 
determination by Scho nrich et al. 1 (|2010h . 

The results for the 2D (Vr) and (V^) velocity maps in 
the Galactic plane for a 4 x 4 kpc box (60 x 60 bins) cen- 
tred on the Sun are presented in Figs. U and [5] for both sets 
of parameters and for both samples (full sample and red 
clump selection). X and Y increase towards the GC and in 

3 Although this does not formally exist in a non-axisymmetric 
potential. 



the direction of Galactic rotation respectively. This places 
the GC either at (8,0) or (7.8,0) depending on the Galac- 
tic parameters used. The solar circle (R = Ro) is drawn as 
a dashed line and the open circle delimits a zone 125 pc 
in radius similar to the Hipparcos sphere. For orientation, 
the location of the local (Orio n) and Perseus arms are indi- 
cated based on the CO map of Engl maier et al. ~l (|2010h . No 
strong kinematic signal is associated with them in (Vr), but 
a strong decrease of (V^) is associated with the Orion arm 
(while little data are available on the Perseus arm). This 
could help place limits on the size of the potential pertur- 
bations associated with these arms, but is beyond the scope 
of this paper. 

Again, the global structure of the velocity field in figs. [4] 
and [5] comparing the top and bottom panels, is similar for 
both samples. Small differences exist on local scales, espe- 
cially around (X,Y) = (0.8,— 0.4), but aside from this region 
the velocity trends are compatible. Isocontours in (Vr) are 
orientated about 50° from the direction to the GC and do 
not follow a particular symmetry axis of our procedure or 
sample. Within the local neighbourhood indicated by the 
open circle, no net motion is detected and our values are 
compatible with previous measurements. The velocity gra- 
dient becomes noticeable only outside the Hipparcos sphere. 

Comparing left and right panels of figs. U and [5] we see 
that using updated Galactic parameters does influence the 
measured mean motions. In particular, in Fig. [4] it dimin- 
ishes but does not eliminate the mean motions. Trying to 
minimise the gradient while keeping the ratio v c o/Ro con- 
stant, leads to unrealistic values for the Galactic parameters 
with Ro < 4 kpc and v c o < 126 km s _1 . 



4 DISCUSSION AND CONCLUSIONS 

In this paper, making use of RAVE line-of-sight velocities 
and distances alone, we report the discovery of a radial veloc- 
ity gradient in the Galactic disc amounting to at least 3 km 
s _1 kpc -1 observed in the direction of the GC for d < 2 kpc. 
Making use of the full sample, we fitted the Oort constants in 
the classical way, and obtained values globally in agreement 
with those of OD03 for old red giants. But as their measure- 
ment was based on proper motions alone, they could not 
constrain the value of the Oort constant K, which we esti- 
mate here to be roughly between 4 and 6 km s _1 kpc -1 (see 
Table 1). 

Then, assuming a zero radial velocity of the Local 
Standard of Rest, we reconstructed the 2D velocity map, 
and found that the Galactic parameters do play a role in 
the measured amplitude of the streaming motion. Adopting 
standard IAU values, which are suspected to underestimate 
Vco/Ro, yields a larger amplitude. On the other hand, us- 
ing more recent values for v c o and Ro reduces the net radial 
motion but does not eliminate it. The latter result is inde- 
pendent of the source of proper motions and distances used. 
We note that the confirmed existence of such a net outwards 
motion for stars, amounting to 10% of the circular speed, 
would raise the question of the validity of determining the 
Galactic parameters b ased on the assum ption of circular mo- 
tion (see for example iReid et al. Il2009h , and the error bars 
on the measured ratio v c o/Ro are likely underestimated. We 
however also note that, in order to establish it more firmly, 
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Table 1. Oort constants A, B, C and K measured using the RAVE sample. The kinematic centre Vv = 0.5 arctan(— C/A), the amplitude 
VA' 2 ±C 2 , and the velocity gradient dV R /dR = K + C are also given. All numbers (apart from %j) v ) are in km s 1 kpc 1 . The first 
two lines give the old values from KT94 and OD03. For the RAVE sample, errors bars are computed using a Monte-Carlo sampling, 
but they do not take into account the possible systematics listed in OD03. The last two lines show changes in the Oort constants if one 
considers a 20% bias in our distances. The two cases are obtained by multiplying the distances of the whole catalogue by 0.8 or 1.2, then 
recomputing the velocities and re-applying the same R, z selection as before. 





Sample 


A 


B 


C 


K 


Vv(deg) 


VA' 2 ± C 2 


K + C 


KT94 


^los, Vb, and in 


14.4 ±1 


-12.0 ±3 


0.6 ± 1 


-0.35 ±0.5 


-1.2 


14.4 


0.25 


OD03 


Vl 


15.9 ±2 


-16.9 ±2 


-9.8 ±2 




15.5 


18.7 




RAVE 


vi os and fi b 


13.6 ±0.5 




-9.6 ±0.5 


5.7 ±0.3 


17.6 


16.6 


-3.9 


RAVE 


v\ os , fi b and iii 


9.2 ±0.5 


-17.4 ±0.5 


-12.7 ±0.5 


4.6 ±0.5 


27.0 


15.7 


-8.1 


d x 0.8 


v\ os and ii h 


15.0 




-9.2 


6.3 


15.8 


17.6 


-2.9 


d x 0.8 


vios, Vb and Vl 


11.1 


-11.8 


-10.5 


5.6 


21.7 


15.3 


-4.9 


d x 1.2 


vi os and ii b 


12.5 




-10.1 


5.3 


19.5 


16.1 


-4.8 


d x 1.2 


vios, Vb and jii 


7.9 


-21.3 


-14.4 


4.1 


30.6 


16.4 


-10.3 



our result should imperatively be confirmed by further line- 
of-sight velocity measurement s at lower latitudes with fut ure 
spectroscopic surveys (see e.g. iMinchev &; Quillen ll2008l ). It 
is indeed very intriguing that a strong 21-cm absorption fea- 
ture along the line-of-sight to the GC is compatible with a 
zero mean velocity relative to the LSR (Radhakrishnan & 
Sarma 1980). 

The cause of such a large-scale streaming motion would 
clearly be a non-axisymmetric component of the Galactic 
potential. It could thus be due to: (i) the Galactic bar; 
(ii) spiral arms; (Hi) the warp; (iv) a triaxial dark mat- 
ter halo; or (v) a combination of some or a ll of these com- 
ponents. Indeed, [Minchev &; Quillenl (|2007h have simulated 
the effect of spiral density waves on the Oort constants and 
have showed that for spiral-only perturbations, the value of 
\C\ is larger for lower stellar velocity dispersions, contrary 
to the findings of OD03. The observed trend in C could, 
however, be explained with the effect of the Galactic bar 
(Min chev etHI I2007T ) , but not its exact magnitude. Actu- 
ally, given that the Milky Way disk contains both spiral 
arms and a central bar, the effect of both these perturbers 
needs to be considered simultaneously. In external galax- 
ies, such as M81, such non-circular motions in gas flows 
are usually associated with spiral arms a nd/or the central 
bar. ISellwood &; Zanmar Sanchez I ([2010) have also placed 
bounds on the ellipticity of the outer dark matter haloes 
of extern al galaxies from the observed non-circular mo tions 
(see also iBinnev |[l978l : iFranx fe de Zeeuw Ill992l : lKT94h . It 
will be of the highest interest to determine which model will 
best reproduce the trend we report in this paper. However, 
it is unlikely that we will be able to discriminate between the 
various models with such a small fraction of the Milky Way 
sampled. More, deeper observations from future astromet- 
ric and spectroscopic surveys will probably be mandatory in 
this respect. 

Observations of non-circular motions of gas in the inner 
Galaxy determine the characteristics of the central region of 
the Milky Way and constrain the parameters of, for exam- 
ple, the Galactic bar and the amount o f dark matter allowed 
inside the solar circle (se e, for instance, Bissant z et al. 1 120031 : 
iFamaev fe Binnev ll^OoBh . It is likely that mapping and un- 
derstanding the non-circular motions in the outer parts of 
the Milky Way will reveal similarly important information 
about the Galactic potential and the mass distribution on 
larger scales. 
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Figure 4. Mean velocity fields (Vr) derived using Eq. [2]and assuming R = 8 kpc and v c o = 220 km s _1 (left panels) and Ro = 7.8 kpc 
and v c o = 247 km s _1 (right panels). Top panels: full sample . Bott om panels: red clump selection. The locations of the nearest spiral 
arms are indicated using the CO map from Englmaier et al. I (2010). The open circle delimitates a sphere 125 pc in radius around the 
Sun. The maps are smoothed over 3 pixels to highlight the mean velocity trends. The maps are 60 x 60 bins in size between -2 and 2 kpc 
along each axis. X increases in the direction of the Galactic centre, Y is positive towards the Galactic rotation. 
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Figure 5. Same as Fig[4]for (V^). 



